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EXPERIMENTAL AND PREDICTED LONGITUDINAL RESPONSE 

CHARACTERISTICS OF A LARGE FLEXIBLE 35° SWEPT- 

WING AIRPLANE AT AN ALTITUDE OF 35,000 FEET 

By Henry A. Cole, Jr., Stuart C. Brown, 
and Euclid C. Holleman 

SUMMARY 


The longitudinal frequency response of a large flexible swept-wing 
airplane as determined from transient flight data, excited by elevator 
pulses, is presented for flight conditions of 0.6 to 0.85 Mach n\mber at 
an altitude of 35^000 feet for a range of center-of -gravity locations. 
These data cover a band of frequencies which include the short-period 
mode and two aeroelastic modes. Response quantities are presented for 
center-of -gravity, wing, and tail locations in order to define the mode 
shapes. 

Predicted transfer functions based on the two-degree-of -freedom 
longitudinal equations of motion, with coefficients modified to include 
flexibility effects, are compared with experimental transfer functions 
evaluated from the measured frequency responses. Comparisons are also 
made with transfer functions predicted for the rigid airplane in order 
to show the effects of flexibility. The results are presented in the 
form of transfer -function coefficients for the response of the center of 
gravity and show good agreement between predicted and experimental values, 
except at conditions where the airplane has a pitch -up tendency. For 
responses near the wing tip, this theory becomes inadequate for frequences 
near and above the wing first-bending mode frequency. 

A method for determining the aerodynamic lift and moment of a flex- 
ible wing through use of aerodynamic and structural influence coefficients 
is presented. 


INTRODUCTION 


The desire to increase the range and speed of large airplanes has 
recently led to configurations with sweptback wings of high aspect ratio, 
thin airfoils, and fuselages of high fineness ratio. All of these factors 
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tend to increase the flexibility of the structrire, and the associated 
aeroelastic effects are becoming of greater importance in problems of 
static and dynamic stability and control. The dynamic effects are espe- 
cially important when the airplane is equipped with automatic control 
because structural modes may introduce system instabilities which would 
not be present in a rigid airplane. It is important to be able to pre- 
dict these aeroelastic effects in order to insure that the basic handling 
qualities and loads requirements are satisfied, and to permit rational 
design of automatic control systems concurrently with the airplane. 

In view of the above problems, the NACA is currently flight testing 
a large flexible 35° swept-wing airplane over its entire operating range 
of Mach numbers and dynamic pressures. The aims of this program are to 
document and analyze the airplane response to control surface motions 
and, through comparisons between measured and predicted response charac- 
teristics, to establish simple but adequate methods of prediction for 
flexible airplanes. The flight data considered in the present report 
consist of measurements of longitudinal response to elevator deflection 
at an altitude of 35,000 feet over a Mach number range of 0.6 to 0.85 
for several center-of -gravity positions. The measured response quanti- 
ties (nomnal acceleration at several stations, pitching velocity at the 
center of gravity, and wing deflections) were suitable for frequency- 
response analysis over a range of frequencies including the short-period 
longitudinal mode and the first two aeroelastic modes of the airplane. 

With regard to test and analytical techniques, one convenient way 
to express response characteristics is in terms of frequency response. 

The longitudinal frequency response may be evaluated directly from flight 
by measirring the response quantities resulting from sinusoidal elevator 
motions of various frequencies. However, a considerable saving in flight 
time is realized if, instead, the airplane is excited by a pulse -type 
elevator motion and the measured transient time histories are transformed 
to frequency-response form by the Fourier integral, a technique described 
in references 1, 2, and 3* In the case of the test airplane, preliminary 
calculations indicated that pilot-applied pulses would provide adequate 
excitation over the range of frequencies of interest and, hence, this 
technique was used. Another application of this method to a relatively 
flexible airplane is reported in reference 4. 

Although these frequency-response data completely define the longi- 
tudinal response of the airplane, operational expressions relating the 
output response to an input distiorbance, known as transfer functions, 
are a more useful form for detailed analysis or for the synthesis of 
automatic control systems. These transfer functions may be approximated 
from the measured frequency response by a curve -fitting procedure such 
as that described in reference 5* Comparable predicted transfer func- 
tions may be derived from the equations of motion including rigid body 
and structural modes. 
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Simplifications in these transfer functions appeared justified for 
the range of flight conditions covered in this report. Consequently, 
theoretical transfer functions were derived from the longitudinal sta- 
bility equations for two degrees of freedom, and first-order effects of 
flexibility were included by modifying the coefficients of these sim- 
plified transfer functions. The coefficients were calculated for both 
the rigid and flexible airplane for comparison with coefficients 
evaluated from flight-test data. 

Data used in this report were obtained from flight tests conducted 
at the High-Speed Flight Research Station of the NACA and the analysis 
and reduction of data were a cooperative effort of HSFRS and Ames Aero- 
nautical Laboratory. 
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NOTATION 

amplitude 

lift coefficient 

pitching -moment coefficient 

differential operator, — 

dt 

imaginary part 
normal-acceleration gain 
pitching-velocity gain 

radius of gyration about principal lateral axis, mean aero- 
dynamic chords 

Mach number 

real part 

wing area, sq ft 

normal acceleration time constant, sec 
pitching velocity time constant, sec 

velocity, ft /sec 
weight, lb 


wing span, ft 
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c.g. 
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wing chord, ft 

center of gravity, percent c 


wing mean aerodynamic chord 


' 1 / 


h/2 


c^dy 


acceleration due to gravity, 32.2 ft/sec^ 
mass, slugs 

normal acceleration, positive downward, gravity units 
dynamic pressure, Ih/sq in. 

longitudinal distance from center of gravity to subscript 
quantity, positive when center of gravity is foivard of 
subscript qtiantity, ft 

spanwise coordinate perpendicular to plane of symmetry, ft 
structural deflection, positive downward, in. 
angle of attack, radians 

elevator control deflection, positive downward, radians, 
except where noted 

change in streamwise angle of attack due to wing distortion, 
radians 

damping ratio, dimensionless 

spanwise coordinate, fraction of wing semispan 

pitch angle, radians 

mass density of air, slugs/cu ft 


m 

PSV' 


sec 


*^/output\ angle of output qiiantity minus phase angle of input 

, quantity 


0) 


frequency, radians/sec 
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undamped natural frequency, radians/sec 


Subscripts 


t 

tail 

wt 

wing tip 

a 

A 

8a 

d 

8 

8(da/dt) 

e 

A 

8e 

e 

8 

8(d©/dt) 

e 

8 

8(de/dt) 
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TEST EQUIPMENT 


The test airplane was a Boeing B-47A with General Electric J-47-GE-23 
turbojets and with wing vortex generators as shown in figure 1. The air- 
plane was fitted with an external nose boom and an optigraph on the top 
of the fuselage which recorded the movement of 100-watt target lights 
located on both wing panels and on the tail (figs. 1 and 2). Elevator 
angle was measured by six NACA resistance-type control-position indicators 
located on the left and right elevators at the root, midsemispan, a n d tip. 
The outputs were recorded on Weston 12-channel and Consolidated l8-channel 
oscillographs. The pitching velocity at the center of gravity was meas- 
ured by a magnetically damped NACA pitch turn meter, the acceleration at 
the center of gravity and tail by NACA air-damped accelerometers, and the 
acceleration at the wing tip by Statham linear accelerometers. 

The locations of the optigraph targets and wing stations used in the 
present study are shown on figure 3. This figure also gives information 
on the distribution of wing weight which is used subsequently in the 
analysis. 
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MEASURED FREQUENCY RESPONSE 
Test Procedure 


The flight-test conditions covered Mach numhers from 0.6 to 0.85 at 
an altitude of 35^000 feet for a range of center-of -gravity locations 
from 12 to 30 percent of the mean aerodynamic chord and a gross weight 
from about 110,000 to 120,000 pounds. Table I lists the flight condi- 
tions used in this report. 

Before each maneuver, the airplane was trimmed at the desired speed 
and altitude. Then the pilot^ applied a quick elevator pulse and held 
the stick fixed until the transient motion damped out. It was found that 
the pilot ■was capable of consistently applying a 4° to 6° elevator pulse 
of approximately 0.5-second duration and returning (nearly) to the trim 
condition. All of these pulses produced about l/4 g acceleration at the 
center of gravity. During these maneuvers records were taken of elevator 
angle, pitching velocity at the center of gravity, normal acceleration 
at the center of gravity, tail, and wing tips, and wing deflections. 
Typical time histories of the elevator control input and the measured 
output quantities are shown on figure 4. 

The' longitudinal moment of inertia was measured by ground oscilla- 
tion tests in which the airplane was supported on knife edges and a 
spring . 


Reduction of Data 


The transient data measured in flight represent the longitudinal 
dynamic response of the airplane to a particular elevator input, whereas 
for detailed analysis it is desirable to know the airplane response to 
an arbitrary input. As shown in references 1, 2, and 3^ it is possible 
to transform the input and output response quantities into frequency- 
response form. This defines the response of the airplane to sinusoidal 
elevator motion of various frequencies. Since the response to any arbi- 
trary input may be obtained by applying an inverse transformation to the 
frequency response, the frequency response defines all of the basic 
characteristics of the airplane response independently of the particular 
input, and in itself provides a means for studying the nature of the 
response and the effect of variables. 

Method . - By means of the Fourier integral the transient function of 
the input or output, f(t), is for zero initial conditions transformed 

^Acknowledgement is given to Aeronautical Research Pilots 
Joseph A. Walker and Stanley P. Butchart. 
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into the frequency function, g(ico), as follows: 


g(iw) = / f(t) e dt 


Separating g(iw) into real and imaginary parts gives: 


R[g(iO))]= / f(t) cos wt dt 


l[g(io))]=- / f(t) sin Lit dt 
^o 


The amplitude and phase angle obtained from the real and imaginary 
parts are: 


A = Vr^ + 


9 = tan"^ i 

The method used for evaluating the above integrals was similar to 
the one described in reference 6, in which ordinates of the time history 
are read at uniform intervals such that the time history is closely 
approximated by a series of parabolic arcs. (All traces were read at 
0.05-second intervals except for the pitching velocity and wing deflec- 
tion which were read at 0.1-second intervals.) The integrals were 
determined by multiplying these ordinates by a set of coefficients for 
each frequency and then summing the products. Since the ordinates of 
the time history were read only until a steady-state value was reached, 
which usually differed from zero, an analytical expression was used for 
the remaining portion of the integral. These operations were calculated 
on card-programmed IBM digital computing equipment. Corrections were 
made in the data for the dynamic response of the instruments . 

Accuracy . - The accuracy with which the frequency response may be 
determined is dependent on the relative magnitudes of the transform of 
the transient quantities and the transform of errors arising from instru- 
ment and reading inaccuracies. These errors in the time history were 
estimated to be of the order of 0.005 inch of film deflection. Therefore, 
for linear calibration curves, the transform of the errors consists of 
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the Foixrier transform of random errors up to about 0.005 inch over the 
length of the transient plus a step error of 0.005 inch in the analytical 
correction at the end of the record. This latter value is by far the 
largest and for practical purposes may be considered to be a measure of 
the maximum expected error. The data in this report were considered to 
be sTifficiently accurate if the amplitude of the Fourier transform of a 
step of 0.005 inch of film deflection was less than 10 percent of the 
amplitude of the total transform. If the transform of the errors is out 
of phase with the transform of the transient, then the error in the 
amplitude will be less than the value above and the phase angle will be 
in error. However, this phase-angle error will be less than 6° if the 
amplitude of the error is less than 10 percent of that of the transient. 
For this reason, accuracy in the amplitude also insures accuracy in the 
phase angles . 

In initial flight tests a nimiber of pilot-applied elevator pulses 
of varying length were recorded and the amplitudes of their transforms 
were compared with the above accuracy criterion as shown on figure 5. 

The error boundary is the amplitude of the transform of a step of 0.6° 
elevator deflection (which corresponds to 0.005 inch of film deflection 
multiplied by 10 ) and represents a boundary of amplitudes below which 
errors of greater than 10 percent can be expected. The frequency range 
was selected to include the first two longitudinal aeroelastic modes 
indicated by ground vibration tests in reference 7, and the natural fre- 
quencies of these modes are marked on this figure. It may be seen that 
the transform of the longer pulse (run A) tends to go to zero near the 
frequency of the wing first-bending mode and falls below the accuracy 
boundary at regular intervals thereafter. Hence, the longer pulse does 
not provide adequate excitation for accurate evaluation of the frequency 
response at frequencies of the aeroelastic modes. The shorter ' pulse 
(run B), on the other hand, provides adequate excitation for frequencies 
from 1 to 20 radians per second which is the range of interest in this 
report. Hence, pulses with a time base of about 0.5 second, as in run 
B, were used to obtain the frequency-response data in this report. 

Frequency-response curves evaluated from transient data in which 
the airplane was excited by pulses of different lengths at the same 
flight condition are compared on figure 6. Also shown on this figure 
are boundaries below which errors of greater than 10 percent would be 
expected. This error boundary represents the value of Iq/s] for which 
the sum of the separate errors of 1©| and |s| is equal to 10 percent. 

It may be seen that the curves agree within 10 percent over the range 
where the accuracy criterion is satisfied. At frequencies below 1 radian 
per second, the frequency response obtained by the shorter pulse becomes 
inaccurate and the curves disagree. Also, at frequencies above 7 radians 
per second, the frequency response obtained by the longer pulse becomes 
inaccurate and it may be seen that the values become erratic. Also 
shown on this figure is a frequency response evaluated from transient 
data excited by a short pulse on another flight in which conditions were 
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slightly diff ©rent • This curve agrees well with the other frequency- 
response curve obtained with a short pulse over the frequency range 
of 1 to 20 radians per second, except that it is displaced upward at' all 
frequencies. As will be shown later, this difference may be accounted 
for by the difference in the time parameter t for the two flights. 

Thus, the data of figure 6 show that with the short-pulse excitation, 
the frequency response can be evaluated with stifficient accuracy over 
the frequency range of 1 to 20 radians per second. 

Since the frequency-response technique is only valid for linear 
systems, frequency responses are questionable when obtained at flight 
conditions in which aerodynamic derivatives are believed to be nonlinear. 
This airplane appears to have a linear response at 0.72 Mach number and 
0*^ lift coefficient since the same frequency response was obtained 
with different degrees of excitation in figure 6. However, a study of 
wind-tunnel data indicated that for flight conditions above 0.8 Mach 
number and for lift coefficients above 0.6 at lower Mach numbers the 
stability derivatives become nonlineajr. In view of this factor caution 
should be exercised in extrapolating the test results obtained at flight 
conditions close to these boundaries (see table I for flight conditions) 
to disturbances greater than those used in the flight tests, that is, 
A nmay = lA 8 at the center of gravity. 

No corrections were made in the frequency responses for effects of 
fuel sloshing. 


Results and Discussion 


Frequency-response curves evaluated by the foregoing method are 
shown on figures 7 and 8 for the test range of Mach numbers at an alti- 
tude of 35,000 feet; quantities measured at the wing tip and tail as 
well as at the center of gravity are included so that a fairly complete 
picture of the response of the airplane is presented. The frequency 
responses are shown for a forward center-of -gravity location in figure 7 
and for a rearward one in figure 8. In the following discussion on 
these frequency-response curves , the various longitudinal modes will be 
identified by comparing the accelerations at the center of gravity, wing 
tip, and tail, and wing deflections. Also, the frequencies at which 
these modes occur will be compared with those obtained by the ground- 
vibration tests of reference 7* 

Short-period mode .- The first and largest peak in the amplitude 
ratio of all of the frequency-response curves is the short-period longi- 
tudinal mode which occurs at frequencies from 1 to 3 radians per second. 
These values agree approximately with frequencies estimated from wind- 
tunnel data. The accelerometer responses at wing tip., tail, and center 
of gravity, which are compared in figure 9, show that the accelerations 
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at the three locations are essentially in phase, and that the amplitudes 
are larger at the ving tip and the tail. The acceleration response at 
these latter stations differs from that at the center of gravity because 
of the contributions of structural deflections and pitching acceleration, 
which will now be discussed. 

The wing -bending deflection response, as determined from optigraph 
records, is presented on figure 10. At the short-period frequency, all 
deflections are similarly phased and increase in amplitude toward the 
wing tip, representing a type of deflection similar to the wing first- 
bending mode. Since for sinusoidal motion acceleration is equal to minus 
the frequency squared times the displacement, the contribution of the 
wing deflection to the amplitude ratio of the acceleration at the wing 
tip is of the order of -1 g per radian. By the same argument, the con- 
tribution of tail deflection to the amplitude ratio of the acceleration 
at the tail was also found to be small. Hence, the contribution of 
deflections to accelerations at the wing tip and tail are small for fre- 
quencies in the neighborhood of the short-period mode. 

The contribution of pitching acceleration to the acceleration ampli- 
tude at the wing tip and tail may be simply calculated by multiplying 
the amplitude of the pitching acceleration (0to) by the distance to the 
center of gravity (fig. 2). This is permissible because the phase angles 
of the acceleration response at the center of gravity and the pitching 
acceleration response (phase angle of O/b equals phase angle of 0/5 
plus 90^) are similar at the frequency of 2 radians per second. The 
increment in acceleration response at the wing tip and tail, then, is 
about 3 and 5-5 g's per radian, respectively, which agrees with the added 
increments in the acceleration peak shown on figure 9* Hence, the 
increase in the acceleration amplitude at the wing and tail is primarily 
caused by the airplane pitching acceleration. 

Comparison of the phase angles on figure 10 with those on figure 9 
shows that the wing deflections are in phase with the accelerations and, 
therefore, the wing tips deflect in proportion to and in the same direc- 
tion as the airplane accelerates (sketch (a)). This sketch illustrates 
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Sketch (a).- Short-period mode. . 

the relationship of the wing bending to the airplane center-of-gravity 
motion. Actually, as shown in the discussion of pitching acceleration, 
the wing tips travel farther in space than the center of gravity because 
of the pitching motion of the airplane. 
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Wing first-bending mode .- The next peak in the frequency response, 
as seen on figures 7> ^3, and 10, occurs at frequencies of from 8 to 9 
radians per second, slightly higher than the wing first-bending frequency 
(6.9 radians per second) in ground -vibration tests. Calculations indi- 
cate that the frequency of this mode is higher in flight than on the 
ground because of the additional spring force contributed by aerodynamic 
forces and the increased freedom of the body in pitch and translation. 

As seen on figures 7(a) through' (c) and 8(a) through (c), the response 
peaks for locations on the fuselage are small for this mode as compared 
to those of the short-period mode, but at the wing tip (figs. 7(d) and 
8(d)) a very high peak occurs. Referring to figure 9, it is noted that 
the wing-tip acceleration undergoes a l80° phase shift at this peak. 

Also, according to figure 10, the amplitudes of the wing deflections 
increase toward the wing tip and are in phase with each other and l80° 
out of phase with the wing-tip acceleration which establishes this as 
the wing first-bending mode. The deflections are also in phase with the 
acceleration at the center of gravity, as was the case for the short- 
period mode. The main characteristic which distinguishes this mode from 
the short-period mode is that the wing deflections are the largest factor 
in the accelerations at the wing tips (sketch (b)), while in the short- 
period mode the body translation is the largest factor. 



Sketch (b).- Wing first-bending mode. 

Other modes .- Several small peaks appear at frequencies from l4 
to 16 radians per second, but these are not well defined because the. 
frequency response is in a region of low accuracy. These peaks are most 
prominent on the wing-tip acceleration responses (figs. 7(d) and 8 (d)). 
The wing-deflection response on figure 10 indicates an upward trend in 
the mid-semispan-deflection amplitude ratio as compared to that of the 
wing tip. This would indicate a mode of the wing second-bending type 
involving body translation and pitch. Some calculations were made on 
the modes of vibration of the B-47 with body translation and pitch 
included; these indicated a mode of the type shown in sketch (c) at I6.7 



Sketch (c).- Wing second -bending mode coupled 
with body translation and pitch 
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radians per second. This agrees qualitatively with the acceleration 
frequency responses (figs. 7(b), (d) and 8(b), (d)), which show that the 
accelerations of the wing tip and the center of gravity tend to become 
more closely phased when the mode becomes prominent at the high Mach 
nxanbers. Another possible mode of vibration in this frequency range 
involves bending of the inboard nacelle supporting structure which, as 
indicated by ground-vihration tests in reference 7, excites considerable 
wing motion at l6.k radians per second. This mode is probably closely 
coupled with the mode mentioned above, which would explain the appear- 
ance of the two closely spaced peaks in this frequency region. 

Effect of Mach n-umber and dynamic pressure . - Since Mach ntimber and 
dynamic -pressure effects cannot be separated when, as in the present 
case, data are available for only one altitude, they will be considered 
together with Mach number arbitrarily selected as the independent vari- 
able. From figures 7 and 8 it may be seen that there is a gradual 
increase in the amplitude of the short-period mode peak for the response 
quantities up to about a Mach number of 0.8l, after which the trend 
reverses, probably because of critical Mach number effects. The general 
level of the response quantities also follows this trend at the higher 
frequencies, although the acceleration at the wing tip seems to be rela- 
tively insensitive to Mach number changes. At these higher Ife,ch numbers 
the frequency response is somewhat dependent on the magnitude of the 
elevator pulse input because of a nonlinear effect mentioned previously. 

The frequency of the amplitude peak of the short-period mode also 
increases up to a Mach number of 0.8l above which it decreases. The 
frequencies of the peaks of the higher modes are relatively constant over 
this range of Mach nvimbers and dynamic pressures, although the same trend 
may be noted. 

Effect of center-of-gravity location .- From figure 11(a) it is seen 
that the effect of moving the center of gravity from 12.6 to 29.7 has 
little effect on the frequency response of pitching velocity, although 
the peaks of the short-period mode occur at lower frequencies as the 
center of gravity is moved back. The effect of center-of ^gravity move- 
ment is more apparent on the acceleration response (fig. 11(b)), which 
shows a definite trend of the short-period mode to higher peaks and 
lower frequencies as the center of gravity is moved back. 

For the rearward center-of-gravity location there are from 3,000 
to 4,000 pounds of additional fuel in the aft main tank (fig. 2) as com- 
pared with the forward center-of-gravity location. However, there 
appears to be no significant change in the high-frequency portion (5 to 
15 radians per second) of these responses. This might be expected, 
since the fuselage vertical-bending mode occurs at a much higher fre- 
quency (29 radians per second) . 
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EXPERIMENTAL AND PREDICTED TRANSFER FUNCTIONS 


Although frequency-response plots completely defined the response 
characteristics over the frequency range of interest, analytical expres- 
sions for the response are more useful in detailed analyses or in the 
syntheses of automatic control systems. Such analytical expressions, 
often called transfer functions, may he evaluated either from the experi- 
mental frequency-response plots or from the predicted equations of motion 
of the airplane . In order to show the relationship between the transfer 
function and the frequency response more clearly, the predicted transfer 
functions will be derived first. Then the method for evaluating experi- 
mental transfer functions will be explained and, finally, the results of 
the experimental and prediction methods will be compared, • 


Predicted Transfer Functions 


Equations of motion .- To define completely the airplane dynamic 
system, the equations of motion must include all of the degrees of free- 
dom, but for practical purposes the degrees of freedom are usually kept 
to a minimum by including only the most significant airplane modes. In 
the present analysis, only the short-period mode is considered since it 
was previously shown that the effects of other modes on the responses of 
pi’incipal interest, those involving motions of the fuselage, were small 
for the test altitude and range of Mach numbers. However, for response 
quantities near the wing tip where at frequencies near and above the 
wing first-bending mode frequency the response is large, these equations 
are, of course, not adequate. 


The equations for longitudinal motion, where changes in forward 
velocity are neglected, were used, with pitching -acceleration (&') terms 
added to account for the distortion of the airplane due to pitching- 
acceleration inertial loads. These equations in operator form aret 


CLa + 2 t^D a + (c^ -2 t)d + C^. D^ 0 + = 0 (l) 


^m^j + a 


Cm^ D+( Cm^ 


i^TKy25 

2V 


)d-] 


® + CniQ 5-0 ( 2 ) 


In order to account for the first-order effects of flexibility in 
the preceding equations, the structural deformation associated with a 
coefficient of the variables a or 0 and their derivatives was assumed 
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to te in phase with the variable, that is, the damping and inertial 
forces due to structural motion were neglected. 

Estimation of stability derivatives .- Stability derivatives were 
derived by available theoiy with the exception of the elevator- 
effectiveness derivatives, which were modified to include additional 
Mach nimiber effects indicated by wind-tunnel tests. The general methods 
applied were those of reference 8 to determine airloads and those of 
reference 9 to account for flexibility effects. Mach nirnber effects 
were included by the Prandtl-Glauert rule as used in reference 8, which 
was indicated by wind-tunnel tests to be not greatly in error up to a 
Mach n\jmber of 0.75* 

In applying the general methods mentioned above, several modifica- 
tions were made to facilitate use of structural test data. The struc- 
tural stiffness of the wing was expressed in the form of influence coef- 
ficients for the front- and rear-spar chordwise locations which are com- 
patible with the measured influence coefficients and with the spanwise 
stations used in determining the span loading in reference 8. These 
wing influence coefficients were obtained from load deflection data pre- 
sented in reference 10. The weight distribution shown on figure 3 was 
separated into equivalent weights as indicated and used with the influence 
coefficients for determining the wing distortion from inertial loadings. 
Fuselage influence coefficients were obtained from the Boeing Airplane 
Company, the airplane manufacturer. Table II lists all the influence 
coefficients, and Appendix A explains how they were used in conjunction 
with aerodynamic influence coefficients to calculate the lift, moment, 
and deflection of the flexible wing and tail due to any initial angle - 
of -attack distribution.- The stability derivatives were calculated by 
determining the initial angle -of -attack distribution due to rigid-body 
motions or to distortion from inertial loads- caused by rigid-body accel- 
erations , and then calculating the resulting lift and moment coefficients 
by the method in Appendix A. The derivatives and the related factors 
which were taken into account are summarized in Appendix B. 

The values of the calculated derivatives are given in table III, 
which lists separately the aerodynamic and inertial contributions. Some 
of the more important derivatives are compared on figure 12 with values 
calculated for the rigid airplane. The Mach number effects are reflected 
by the rigid-airplane derivatives. The flexible -airplane derivatives 
include both Mach number and flexibility effects and show that flexibil- 
ity tends to have the opposite effect of Mach mamber. Of all the deriva- 
tives the largest variation occurs in 0^ , but, in general, aeroelastic 
effects are not exceptionally large because the range of dynamic pres- 
sures in these tests is not large. 
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Evaluation of theoretical transfer functions .- From equations (l) 
and (2), the following transfer function for pitching velocity may be 
obtained: 


0 

6 


1 + 




u)n 


Wn-= 


(3) 


and, from the relation n 


- ^ (e - a) 
g ' ' 


n 

6 


Kn(l + Tn^D + Tn^D^) 


1 + ^ D + ^ 


“n 




(4) 


The acceleration responses at the wing tip and tail are obtained by 
adding the contributions of acceleration at the center of gravity, the 
pitching acceleration, and the structural deflection, so that 


and 


where 



(5) 


(6) 


(7) 
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t = 

J r\ 


'n 


- S-) " "-g (% - 
-% * C„^(cLe -2 t) 


(9) 


J'e = 


-Ci^ Cm^ •^ 

-Ci^ + ^a(^L0 -2 t) 


(.10) 


= 


^a S + 2'^) 

-Cl^ + Cl^ 0^5 


( 11 ) 


K.=-|K; 


( 12 ) 


-Cig CiUq^ + CLq^ Cmg 
-0];^ C]ji^ + Cl^ Cu]g 


(13) 


(li^) 


Since the transfer-function coefficients are complex combinations 
of the stability derivatives and mass parameters, certain terms have 
been grouped together so that the coefficients are analogous to the famil- 
iar equations for rigid airplanes. The derivatives and Cl^ in the 

terms - 2t and + 2t represent primarily the change in effective 

mass of the airplane resulting from lift due to structural distortions. 

As seen in table III, the aerodynamic contributions to these derivatives 
are relatively small. The derivative in the term (4TKy^c/2V) - 0®- 



NACA RM A5lt-H09 


17 


represents the change in effective moment of inertia due to pitching 
moments arising from structural distortions induced by rotational iner- 
tial loads. 

Using the stability derivatives given in table III and figure 12, 
the transfer-function coefficients were calculated for both the flexible 
and rigid airplane and are shown on figure 13. . The coefficient is 

not shown because its value is insignificant. 

Effect of flexibility on transfer-function coefficients .- The 
transfer-function coefficients for the rigid and flexible airplane will 
now be compared to show for the present tests the significance of addi- 
tional terms included to account for flexibility. 

The over-all effect of flexibility on tOn shown on figure 13 
is to reduce the natural frequency by a maximum of about 10 percent. 

This is principally caused by the reduction in Cm^ due to flexibility 
of the fuselage. This effect is partially compensated by a sig- 
nificant increase in the term, which provides an effective 

spring force resulting from the pitching moments arising principally 
from wing deflection due to normal acceleration. The term C’mv in the 
denominator tends to increase the frequency by about 5 percent. 

The effect of flexibility on damping ratio, is small as shown on 
figure 13. The damping forces are reduced by flexibility, but the spring 
and inertial forces are also reduced to such an extent that the over-all 
effect of flexibility on damping ratio is insignificant. 

The pitching-velocity gain. Kg, of the flexible airplane is somewhat 
higher than that of the rigid airplane (fig. 13). This change is prin- 
cipally due to the reduction in magnitude of Chiq^ in the denominator 
which also is in part compensated by the increase in Because of 

this large contribution of • for the flexible airplane, increments 

C7 

in gain for a 10-percent increase in t given on figure 13 show that 
the flexible airplane has about a 15-percent smaller change in gain due 
to changes in gross weight and altitude. 

Flexibility has little effect on as shown on figure 13 . 

Although the reduction in due to flexibility tends to increase 

this is compensated for the most part by the effective reduction in mass 
in the numerator term. 

The discussion of applies to since these factors are related 
by a simple constant. 

The term Tn^ is exceptionally small and is not significant over 
the range of frequencies considered in this report- 
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Flexibility tends to increase Tn^ (fig. 13). Although this term 
is small, it does have a noticeable effect at frequencies greater than 
that of the short-period mode. 

Effect of approximate equations on transfer-function coefficients . - 
Approximate equations are often used for calculating the transfer- 
function coefficients. These equations, which are obtained by simplify- 
ing equations 8, 9, 10, and 11, are as follows: 



(15) 


/ -PSc 
V 32mKy2 J 
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(16) 
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2t 


(17) 


Te 


2t 

CLa 


(18) 


Transfer-function coefficients for the rigid and flexible cases 
were calculated using both the complete and the approximate equations 
and the results are shown on figure l4. 

Coefficients evaluated by the approximate equations show reasonably 
good agreement with those from the complete equations for the rigid air- 
plane, but large discrepancies are apparent in the case of cog and 
for the flexible airplane. This is principally due to the fact that^the 
approximate equations reflect only the large reduction in magnitude of 
*^mg lue to flexibility while in the complete equations, the effect of 
this reduction in magnitude of Ciq^ is partially compensated by a sub- 
stantial increase in the term -CL^Cjn^. Hence, it may be seen that 0^. 

should not be neglected in calculating the frequency and gain of a flex- 
ible airplane. 
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Experimental Transfer Functions 


Experimental transfer functions may be evaluated by applying various 
curve -fitting procedures to the time histories or to the frequency 
responses. Curve fitting of a time history may be accomplished by a 
least-squares fitting method such as described in reference 11. Curve 
fitting of the frequency response^ as in the following analysis, may be 
done by use of a special set of templates described in reference 5 - 
This method may be explained briefly as follows. It can be shown that 
an expression for frequency response is obtained by replacing the dif- 
ferential operator, D, in the transfer function by the frequency vari- 
able, io). The resulting complex number can be factored into first- and 
second-order terms expressed in polar (amplitude and phase angle) form. 
The method of reference 5 involves fitting the measured frequency 
response by graphical addition of templates selected from a set of curves 
which represent a wide range of first- and second-order factors. 

Evaluation of transfer-function coefficients .- As shown previously, 
the accuracy of the frequency response at low frequencies is question- 
able. For this reason, template fitting of the frequency response was 
only considered to be valid for frequencies above 1 radian per second. 
This limitation made it difficult to fit the frequency response by a 
unique combination of templates for all of the nxxmerator and denominator 
terms which are involved simultaneously in the frequency-response expres- 
sion. Because of this, the natural frequency and damping ratio were 
determined by a least-squares curve-fitting method (similar to that of 
reference 11 ) of the pitching-velocity time histoi^’" o\^er the portion of 
the record in which rate of change of elevator position could be neg- 
lected. The appropriate templates for these values of damping ratio 
and natural frequency were then fitted to the center-of -gravity accel- 
eration frequency-response curves to obtain the acceleration gain, Kn> 
and the time constant, The gain was then determined from K^- 

By use of this value of and the previously determined values of ^ 

and (jj, the pitching-velocity frequency response was fitted to determine 
T 0 . A typical template fit is shown on figures 7(a) and (b) for the 
0.63 Mach niimber curve. 

Transfer-function coefficients evaluated from experimental data for 
the forward and rearward center-of -gravity locations are plotted on fig- • 
ure 15 . These values were corrected to a common altitude and gross 
weight in accordance with predicted variations, but in general these 
corrections were small. 


Comparison of Experimental and Predicted Transfer Functions 


Predicted transfer-function coefficients for the rigid and flexible 
airplane are also shown on figure 15 for comparison with the experimental 
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values. The moment of inertia used in the predicted coefficients was 
determined from ground-oscillation tests for the basic airplane, with 
corrections made for condition of the fuel tanks. The stability deriva- 
tives of table III were employed. 

Effect of Mach number .- Referring, to figure 15, it may be seen that 
the experimental values of natural frequency show good agreement with 
the predicted values for the flexible airplane and fall about 10 percent 
below those predicted for the rigid airplane. An exception is noted at 
a Mach number of 0.84 where the measured frequency decreases sharply, 
probably due to the large decrease in as the airplane approaches 

the pitch -up. Also on this figure, it is seen that the scatter in the 
measured damping-ratio points prohibits definite confirmation of the 
predicted small variations with Ifech number, center-of-gravity location, 
and flexibility. However, the general level of the values shows good 
agreement with theory. 

On figure 15 good agreement is indicated for the pitching-velocity 
gain except at Mach numbers of 0.8 and higher, where large unpredicted 
increases in the measured gain are apparent. This again is probably the 
result of the sudden decrease in the static margin as the critical Mach 
number is approached as previously mentioned. On this figure, good 
agreement also is indicated for the pitching-velocity time constant, T^, 
although the experimental values tend to fall consistently a small amomt 
below the predicted ones. 

Effect of center-of-gravity location . - On figure l6, the variation 
of natural frequency and damping ratio with center-of-gravity location 
Is shown for several Mach numbers. The measured values of natural fre- 
quency show good agreement with the frequencies predicted for the flex- 
ible airplane over the test range of center-of-gravity locations. The 
frequencies predicted for the rigid airplane are consistently higher 
than those for the flexible airplane, but show about the same variation 
with center-of-gravity movement as for the flexible airplane. The 
measured and predicted damping ratios show fairly small changes with 
center-of-gravity location. 

In summary, the longitudinal response of the airplane may be ade- 
quately predicted by the theory used herein to account for flexibility 
up to a Mach niomber of 0.8. The largest errors occur in the natiiral 
frequency if flexibility is not taken into account, but in general all 
of the transfer-function coefficients show only small changes over the 
airplane range of dynamic pressures at 35>000 feet. 

Since the transfer-function coefficients are relatively insensitive 
to flexibility effects at these flight conditions, a more sensitive 
transfer function, the wing-tip deflection response, was calculated 
from equation (7) in order to provide a more accurate check on the aero- 
elastic calculations. The method of calculating deflections is shown 
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in Appendix A. The comparison "between experimental and predicted values 
in frequency-response form is shown in figure 17. It may "be seen that 
the amplitude ratios agree very well and that. the phase angles of the 
experimental response tend to lag the predicted ones 20° to 30 ° up to 
frequencies of 5 radians per second. From this, it would appear that 
• the theoretical calculations of the contribution of the wing to the sta- 
bility derivatives are accurate except for a small phase lag which is 
probably due to the inertial and damping forces due to wing bending 
motion which were neglected in the theory. 


CONCLUSIONS 


The evaluation of the dynamic response of a large flexible airplane 
to elevator pulses at 35^000 feet over a Mach number range of 0.6 to 0.85 
and the comparisons with predicted response have led to the following 
conclusions: 

1. The pulse technique provides dynamic -response data that are 
siifficiently accurate for evaluation of frequency response from 1 to 20 
radians per second, a frequency range which includes the first two aero- 
elastic modes. 

2. Except at conditions where the airplane has a pitch -up tendency, 
good agreement between predicted and experimental transfer functions is 
obtained for frequencies near that of the short-period mode, using pre- 
dictions based on the usual two-degree-of -freedom form of transfer func- 
tion in which the coefficients were modified to include zero -frequency 
aeroelastic effects . 

3. For frequencies near and above that of the wing first-bending 
mode, the modified two -degree -of -freedom theory is inadequate. However, 
the response amplitudes at the center of gravity and tail which are asso- 
ciated with airplane structiaral modes are small compared to the response 
amplitude of the short-period mode and could be neglected for this range 
of flight conditions. On the other hand, the response amplitudes near 
the wing tip become relatively high at these frequencies and could not 

be neglected. 

4. For the flight-test conditions, the principal effects of air- 
plane flexibility are to decrease the natural frequency and increase the 
steady maneuvering acceleration per unit elevator deflection, trends 
which may be attributed principally to fuselage and tail bending and the 
associated loss in the angle -of -attack moment derivative, Cm^j. This 
effect would be more significant if it were not compensated partially by 
the effect of the increase in the rate -of -pitch moment derivative. Cm a, 
due to normal acceleration. 
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5 . The method of employing aerodynamic and, structviral influence 
coefficients in aeroelastic calculations which was used in this investi- 
gation proved to be advantageous in simplifying spanwise loading calcu- 
lations and in applying static-test load data. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 9 , 195^ 
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APPENDIX A 

A METHOD FOR DETERMINING THE AERODYNAMIC LIFT AND MOMENT OF 
A FLEXIBLE WING THROUGH USE OF INFLUENCE COEFFICIENTS 


The approach used is generally that presented in reference 9. How- 
ever, for the calculation of wing deformations, the use of structural 
influence coefficients in conjunction with aerodynamic influence coef- 
ficients will be introduced in place of the wing bending and torsional 
stiffness distributions El and GJ, and the distributed aerodynamic load- 
ings used in reference 9* 

The loading on a flexible wing may be separated into two parts: 
that of the rigid wing and that produced by the wing deflection. The 
discussion that follows will be concerned with the determination of the 
loading produced by a wing deflection corresponding to a given rigid- 
wing loading. The change in streamwise angle of attack due to wing 
deformation may be expressed in a power series form as 


e(3) = eo(ri)q + €^(Ti)q2 + €2(Tl)q3 + • • • (Al) 

where 

€Q(q) is the angle -of -attack distribution due to rigid-wing loadirxg 

ei(n) is the angle -of -attack distribution due to loading obtained 
from eo(T)) 

is the angle -of -attack distribution due to loading obtained 
from 

The incremental angle -of -attack distributions €o(ti), may be con- 

verted into incremental lifts to form the series 


ACl = q + Clj q^ + Cl^ q=3 + . . . 


As shown in reference 9^ this series can be represented very nearly by 
the equation 

% 


ACl » 


1 + kq 


(A3) 



2k 
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where 

• k = S a _ — ik , etc. 

The total lift coefficient for the flexible wing may be written as 



(A4) 


where Clj^ is the lift coefficient for the rigid wing. Similar equations 
may be written for the aerodynamic moment coefficient. 

Thus it remains to determine Cl^/Clpj and k through use of the 
influence coefficients. The aerodynamic influence coefficients were 
obtained in the form of the loading coefficient Gyn at a station n 
due to a unit angle of attack at station Vf the angle of attack at the 
other stations being zero. The loading coefficients were obtained by 
the method of reference 8. By use of this method, aerodynamic influence 
coefficients at four spanwise stations could be found which, with proper 
care, were sufficient to provide desired accuracy in determining wing 
deflections. These coefficients, Gyn> were obtained by assuming a unit 
angle of attack at one control point, and zero angle of attack at the 
remaining three, and then solving the set of four simultaneous equations 
consisting of the ayn coefficients obtained from reference 8. For 
instance, calculated from the following simultaneous 

equations 


1 - aiiG^i + a-isGis + ^ 13^13 + ®’14^14 


0 + ^ 22*^12 ®' 23*^13 ®’ 24*^14 
° = ^ 31^11 + ^ 32^12 + ®' 33^13 ®' 34^14 

0 ^ a4iGii + + a-43‘^13 + ^44®14 


^In reference 9f the constant k was determined by the ratios of 
the incremental deflections 

k = = = 5 , etc . 


However, for the present method, it will be more convenient to deal 
directly with the incremental lifts. 
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where subscripts i, 2, 3^ and 4 refer to the semispan stations t| =0.924, 
O.TO 7, 0.383# and 0. The resulting lift coefficient for each xmit angle 
of attack can then he calculated using the equation 


Clv = ^(gv4 + i.848Gv3 + 1.4i4Gv 2 + o.765Gvi) 


(A5) 


The experimentally determined structural influence coefficients (ref. lO) 
were measured in the form of deflections in inches at front and rear 
spars due to 1,000-pound loads at front and rear spars at a ntunher of 
spanwise stations. These were cross-plotted to obtain the influence 
coefficients at the spanwise stations, shown in table II. For use with 
the aerodynamic influence coefficients, these coefficients were foirther 
reduced into the form of a change in streamwise angle of attack in 
radians, due to a 1,000-pound load at the quarter-chord position, Smn . 

The influence deflections due to a unit angle of attack at a con- 
trol station can now be calculated through use of the aerodynamic and 
structural influence coefficients. In order to calculate the deflec- 
tions with siifficient accviracy, an integration formula given in refer- 
ence 8 was used. 


I 

dq = f(q„) sin <Pn 

n=i 

where cp^ = and f(qjj) is the value of f(q) at q = cos For 

the particular n\jmber of control stations used in the present case, and 
since fCq^) = 0, and using only the interval 0 < q < 1, the integration 
formula becomes 



where 


f f(Tl) dq 


3 

Y 

n=i 


(a6) 


11 = 0.1502 

12 = 0.2776 

13 = 0.3628 
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Thus, the angular deflection at station m due to loading due to angle 
of attack at station V, with q = 1 psi is calculated from the equation 


(^o) 


b^(l2)^ — 


1000 


^n^mn^n 


(A7) 


n=i 


with 


m = 1, 2, 3 

V = 1, 2, 3, ^ 

The resulting loading due to (en)., is calculated from the equation 

Vm 


(Ga) 


vn 


3 

z 

m=i 


(^o)vm^^ 


(A8) 


from which the lift coefficient may he obtained 



(Ga)v4+ i*8^(Ga)v3 + o.765(Ga)^^ 


(A9) 


The influence lift coefficient 


loading coefficient (Ga) 

/ \ ^ Vn 

equation for (Gx^j is 




due to deflections due to the 


can he calculated in a similar manner. The 


(cLb)^ =^[(Gb)v4 + 0-765(Gb)^^ 


(AlO) 


Thus, for a given angle -of -attack distribution, av, the lift result- 
ing from the initial angle of attack and from the first and second twist 
distributions is calculated from the following equations 



^ CL^ay 
V=1 


(All) 
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^La = Y. °^V (A12) 

V=1 

4 

*^Lb = Z (^^b)v «V (A13) 

V=1 

The total lift for the flexible wing may then he expressed as 


where 



(A14) 


A similar procedure is used to obtain the aerodynamic moment for a 
flexible wing. 


Wing deflections may also be determined by means of wing-deflection 
influence coefficients in a manner similar to that which was used for 
determining the change in streamwise angle of attack. The deflection at 
station m due to the loading due to angle of attack at station v, 
with q = 1 psi is calculated from the equation 


with 


(^o)vin 


b^(12)^ Y 

1000 Zj 

n=i 


^vn^mn^n 


m = 1, 2, 3 
V = 1, 2, 3, ^ 


(A15) 


where Zmn is the deflection in inches at station m due to a 
1,000-pound load at the quarter chord at spanwise station n. 
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For an arbitrary angle -of -attack distribution, the deflection at 
station m due to the initial load is 


4 



v=i 


(Al6) 


The total deflection for the flexible wing may be developed in a 
manner similar to that used for equation (Al4). Then 


Zm 


q 

l+kq 


(zo)m 


(AIT) 


in which the k from equation (Al4) may be used with sxifficient 
accuracy. 
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APPENDIX B 

EVALUATION OF STABILITY DERIVATIVES 
^La' 


Wing lift-curve slope .- This was determined theoretically from 
reference 8 and Appendix A using the section lift-c\irve slope of 5. 71 as 
determined from wind-tunnel data supplied hy the Boeing Airplane Company. 
This derivative should not be confused with the one in reference 12, 
which includes inertial effects. The present derivative is the one 
which would be evaluated experimentally in a wind tunnel with a flexible 
model. 

Tail lift-curve slope .- This was detennined in a manner similar to 
that of the wing, except that fuselage bending was included in the flexi- 
bility of the tail. It was found that the principal reduction in the 
tail lift-curve slope was caused by fuselage bending. Other factors 
included were the rate of change of downwash (ref. 13 ) and a tail effi- 
ciency factor of 0 . 95 . 

Body and nacelles lift-curve slope .- This was determined from wind- 
tunnel data supplied by the Boeing Airplane Company. 


Ct . 
•Mi 



Lag in wing downwash . - Only the tail contribution was considered 
and was determined in a manner similar to that described in reference l4, 
including aeroelastic effects determined by the method of Appendix A. 
Although this term has small effect from the standpoint of lift, it is 
important in the calculation of Cm-. 


Normal acceleration .- Since normal acceleration is related to d 

in the equation n = - ^(0-d), effects of structural deflections due to 

normal acceleration were included in derivatives in a and 9 by the 
method of Appendix A. The principal contribution to this derivative is 
from the lift of the wing. Contribution of the tail varies according to 
the amount of fuel in the rear main tank and is from 10 to 20 percent of 
the total. 
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Lift arising from angle -of -attack distribution due to pitching 
velocity (primarily a tail contribution) was determined using the method 
of Appendix A. 

Normal acceleration .- Same as normal acceleration part of C^., but 
of opposite sign. “ 




Lift resulting from angle -of -attack distribution caused by struc- 
tural distortion due to rotational inertial loads was determined by the 
method of Appendix A. Primarily a wing contribution, but the total 
effect is small. . 




Rigid-airplane value was obtained from low- speed wind-tunnel data 
supplied by the Boeing Airplane Company and was assxuned to be ’constant 
with Mach number. This is justified because the theoretical increase 
according to the Prandtl-Glauert rule is usually compensated by the 
pressure losses at the elevator hinge point (ref. 15). Aeroelastic 
effects of lift and moment on body bending and of lift on stabilizer 
distortion were included. Stabilizer distortion due to elevator pitch- 
ing moment and elevator distortion were neglected. 


Cm« 


This derivative was determined by multiplying the preceding C];^ 
derivatives by the distance in mean aerodynamic chord lengths from their 
theoretical centers of pressure to the center of gravity. The value 
of Cin^j for body and nacelles was obtained by subtracting the theoreti- 
cal Cnict alone from wind-tunnel values of Cm^j for wing, body, 

and nacelles. It should be noted that this also includes the change 
in Cm^ due to the difference between theoretical and experimental 
of wing alone, a difference which is principally due to a somewhat 
higher loading near the root for the experimental than for the theoreti- 
cal case. Since this additional loading occurs near the wing root, it 
does not affect the aeroelastic calculations and, hence, is appropriately 
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added to in the form of a correction for hody pitching moment. 




Cm^^ 




Cmc, 


These momen't derivatives were determined *by multiplying the corre- 
sponding lift derivatives by their moment arms in a manner similar to 
that described for 
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TABLE I.- FLIGHT-TEST CONDITIONS 


Flight 

number 

Run 

number 

Altitude 

Mach 

no. 


• w 

c.g. 

CM 

T 


16 

36,170 

0.68 

0.56 

121,400 

20.9 

2.4 

5.8 


l8a 

35,500 

.72 

.48 

120,900 

20.9 

2.4 

5.2 


l8b 

35,500 

.72 

.48 

120,900 

20.9 

2.4 

5.2 


20 

36,260 

.78 

.41 

120,000 

21.1 

2.4 

4.9 

2 

21 

34,860 

.82 

.38 

119,800 

21.1 

2.4 

4.4 


22 

34,860 

.84 

.33 

119,300 

21.1 

2.4 

4.3 


23 

34,100 

.85 

.28 

118,200 

21.8 

2.4 

4.1 


25 

32,230 

.63 

.49 

113,600 

21.0 

2.5 

4.9 


26 

35,950 

.62 

.61 

113,200 

21.0 

2.5 

5.8 


15 

35,070 

.63 

.63 

116,100 

12.7 

2.5 

5.6 


l 6 

35,090 

.66 

.55 

115,900 

12.6 

2.5 

5.3 


17 

35,360 

.71 

.46 

115,900 

12.6 

2.5 

5.0 

3 

18 

35,060 

.74 

.45 

115,900 

12.5 

2.5 

4.7 


19 

34,980 

.77 

.38. 

115,900 

12.4 

2.5 

4.5 


20 

34,960 

.81 

.37 

115,700 

12.5 

2.5 

4.4 


21 

35,150 

.82 

.36 

115,600 

12.5 

2.5 

4.3 


22 

35,600 

.84 

.36 

115,500 

12.3 

2.5 

4.3 


1 

35,840 

.79 

.40 

119,100 

24.9 

2.4 

4.8 


3 

34,170 

.73 

.43 

118,200 

24.8 

2.4 

4.6 

5 

4 

34,070 

.69 

.48 

117,800 

24.6 

2.4 

4.9 


5 

34,220 

.65 

.54 

117,000 

24.6 

2.4 

5.3 


6 

35,100 

.60 

.66 

116,600 

24.6 

2.4 

5.9 


15 

35,950 

.72 

.45 

111,100 

21.6 

2.4 

4.9 


19 

34,860 

.73 

.43 

116,200 

29.7 

2.4 

4.8 

Q 

22 

34,590 

.60 

.60 

115,300 

29.5 

2.4 

5.8 

y 

20 

34,780 

.70 

.43 

116,200 

29.7 

2.4 

5.0 


l 6 

34,960 

.81 

.37 

116,800 

29.8 

2.4 

4.4 
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TABLE II.- STRUCTURAL INFLUENCE COEFFICIENTS 


Wing 


1000 -pound load at station^ 

3 F^ 

3 R^ 

2 F 

2 R 

IF 

IR 

5 F ■ 

5 R 

Inches 
deflec- 
tion at 
station^ 

3 F 

3 R 

2 F 

2 R 

IF 

IR 

0.0659 

.0705 

.1769 

.l 8 o 4 

.2302 

.2330 

0.0580 

.0931 

.2052 

.2318 

.2730 

.2940 

0.1765 

.2116 

.7419 

.7881 

1.19^3 

1.2430 

0.1663 

.2398 

.7795 

.9005 

1.2934 

1.4057 

0.2447 

.3057 

1.2134 

1.3309 

2.3486 

2.4598 

0.2307 

.3339 

1.2490 

1.4300 

2.3840 

2.5920 

0.2033 

.2486 

.9285 

.9995 

1.6374 

1.7086 

0.1936 

.2784 

.9590 

1.1070 

1.7184 

1.8677 

Fuselage 

Load 

Change in stabilizer angle 

Ig Normal acceleration 

1 Radian/sec^ pitching acceleration. 

1 , 000 -pound tail load (down) .... 
1,000 inch -pounds moment 

applied at stabilizer 

-0.45° 

p0.521° (c.g. = 17 percent), 
[-0.506° (c.g. =25 percent) 
+ 3 . 42 ° 

0.000207° 


^Station locations are shown on figure 3 * 
tion. R denotes rear spar location. 


F denotes front spar loca- 
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TABLE III.- PREDICTED LONGITUDINAL STABILITY DERIVATIVES FOR THE FLEXIBLE 
AIRPLANE AT AN ALTITUDE OF 35^000 FEET; W = 100,000 POUNDS 


Quantity 

c.g. 


Mach 

number 



0 

0.5 

0.6 

0.7 

0.8 



5.27 

5.31 

5.36 

5.45 

5.6 

Aerodynamic 

CLa 

— 

T.l/V 

.015 

.013 

.011 

.010 

Aerodynamic 

0.12 

36/v 

.069 

.056 

.046 

.038 


.25 

27.2/V 

.050 

.o4o 

.032 

.025 

Inertial 

.12 

.026 1 
g 

.348 

.401 

.444 

.473 

Ct. and -Ct . 

cx 

.25 

V 

.027 - 
g 

.373 

.431 

.480 

.512 


.12 

-.03h 

-.031 

-.030 

-.029 

-.028 

.25 

-.030 

-.027 

-.026 

-.025 

-.024 


.12 

-1.68 

-1.81 

-1.82 

-1.82 

-1.81 

.25 

-1.01 . 

-l.l4 

-1.15 

-1.11 

-1.08 

Aerodynamic 

.12 

-27/v 

-.053 

-.047 

-.o4o 

-.036 


.25 

-26/v 

-.050 

-.045 

-.039 

-.035 

Aerodynamic 

.12 

-lk6/v 

-.302 

-.256 

-.224 

-.204 

■Cm^ 

.25 

-138/v 

-.288 

-.244 

-.213 

-.192 

Inertial 

.12 

V 

.024 - 

.337 

.398 

.451 

.489 

and 

.25 

.024 - 
g 

.338 

.400 

.457 

.499 

Cm^- 

.12 

.031 

.030 

.029 

.028 

.028 

.25 

.030 

.029 

.028 

.027 

.027 

% 

— 

.327 

.312 

.306 

.299 

.292 

Cms 

.12 

.25 

-1.21 

-1.17 

-1.16 

-1.12 

-l.l4 

-1.09 

-1.11 

-1.07 

-1.08 
-1.05 . 

T 

— 

2940/V 

6.05 

5.03 

4.31 

3.77 , 
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Figure 1.- Photograph of the test airplane. 
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Figure 3*~ Wing stations and weight distribution used in analysis 
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Figure 4 .- Typical time histories of input and output response quantities. 
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Figure 6.- Comparison of pitching -velocity frequency responses obtained 
from a long- and a short-pulse elevator input at the same flight con- 
dition and from a short-pulse input at a slightly different flight 
condition; M = 0.72. 
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(a) Pitching velocity at the center of gravity. 


Figure 7 »- Frequency response for the center of gravity at approximately 

12.5 percent c. 
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Frequency, oj, radians/sec 


("b) Acceleration at the center of gravity. 


Figure 7«- Continued. 
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(c) Acceleration at the tail. 



Figure J.- Continued. 
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(d) Acceleration at the wing tip. 
Figure 7*- Concluded. 
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Frequency, w, radians/sec 

(a) Pitching velocity at the center of gravity. 
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Frequency, w, radians/sec 
(c) Acceleration at the tail. 
Figure 8.- Continued. 
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Frequ 


Figure 10.- Wing-deflection 



NACA RM A5te09 



(a) Pitching velocity at the center of gravity. 

Figiire 11.- Comparison of frequency responses for various center-of- 
gravity locations for approximately the same flight condition. 
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Figure 11.- Concluded. 
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Figure 12.- Predicted effect of flexibility on some of the longitudinal 
stability derivatives at an altitude of 35>000 feet. 
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Figure 15.- Concluded., 
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Figiire l6-- Comparison of predicted and experimental variation with 
center-of -gravity position of natural frequency and damping ratio of 
the short-period mode; W = 115^000 pounds. 
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Figure 1T-- Comparison of experimental and predicted wing-tip-deflection 
frequency response for Flight 3^ Run 17 • 
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